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ABSTRACT 

We present precise phase-connected pulse timing solutions for 16 7-ray-selected pulsars re- 
cently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope 
plus one very faint radio pulsar (PSR J1124— 5916) that is more effectively timed with the LAT. 
We describe the analysis techniques including a maximum likelihood method for determining 
pulse times of arrival from unbinned photon data. A major result of this work is improved posi- 
tion determinations, which are crucial for multi- wavelength follow up. For most of the pulsars, 
we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray 
counterpart associations. We report glitches measured in PSRs J0007+7303, J1124— 5916, and 
J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that re- 
veals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able 
to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that 
can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw 
component. Another Chandra ACIS image of PSR J1732— 3131 reveals a faint X-ray point source 
at a location consistent with the timing position of the pulsar. Finally, we present a compilation 
of new and archival searches for radio pulsations from each of the 7-ray-selected pulsars as well as 
a new Parkes radio observation of PSR J1124— 5916 to establish the 7-ray to radio phase offset. 

Subject headings: Gamma rays: stars, pulsars: general, Radio continuum: stars, X-rays: stars 
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1. Introduction 

Pulsar timing involves making precise measure- 
ments of pulse times of arrival (TOAs) at an obser- 
vatory (or spacecraft) and then fitting the param- 
eters of a 'timing model' to those measurements. 
This powerful technique enables extremely high 
precision measurements that probe numerous top- 
ics in fundamental physics and astrophysics. This 
is due to the ability to construct a coherent tim- 
ing model that accounts for every rotation of the 
neutron star over periods of years. Precise tim- 
ing measurements on radio pulsars have yielded 
many fundamental advances including the first in- 
direct detection of energy loss due to gravitational 



timing (Camilo et al. 2002). The timing mod- 



radiation (Taylor et al. 1979) and confirmation 



of many effects predicted by General Relativity 
( |Stairs|[2003l |Kramer fc We^fiW9 1. 

Until recently, pulsar timing was only practical 
in the radio and, in some cases, soft X-ray bands 



( Jackson fc Halpern|2005 Livingstone et al.f2 009, 
for example). For radio and X-ray quiet/faint 
pulsars discovered with the Large Area Telescope 
(LAT) on Fermi, the only option is to time them 
directly using the 7-ray data. Earlier instruments, 
such as EGRET on the Compton Gamma-Ray Ob- 
servatory, required very long exposures to even de- 
tect a handful of 7-ray pulsars and it only observed 
them occasionally, typically during a few 2-week 
viewing periods spread over the 9-year mission. 
With the LAT, we have a vastly more powerful 
instrument for long-term pulsar studies. First, its 
effective area (~ 8000 cm 2 at 1 GeV), energy cov- 
erage (20 MeV to > 300 GeV) and point spread 
function (~ 0.8° at 1 GeV) are greatly improved, 
providing a large increase in instantaneous sensi- 
tivity over EGRET (|Atwood et ~aLl|2009|). Sec- 



ond, because Fermi operates in a continuous all- 
sky survey mode with a very large field of view 
(~ 2.4 sr), it accumulates data on all pulsars in 
the sky roughly uniformly at all times. This al- 
lows long evenly-sampled timing observations of 
all pulsars detectable with Fermi. 

In this paper, we describe the techniques devel- 
oped for precise timing of pulsars using the 7-ray 
photon data provided by the LAT. We then apply 
this method to the first 16 7-ray-selected pulsars 



discovered in blind searches of LAT data (Abdo 
et al.|2009a I plus one additional radio pulsar (PSR 
J1124— 5916), which is too faint for routine radio 



els presented here are updated versions of those 
used for these pulsars in the First LAT Catalog 



of Gamma-ray Pulsars (Abdo et al. 2010d), and 



this paper documents the methods used to cre- 
ate those models. In the case of the bright Vela 
and Geminga pulsars, these methods were used 
to provide high-precision timing models used for 
phase-resolved analysis ( Abdo et al.||2010e|b |. 

Timing observations provide a wealth of impor- 
tant information critical to the understanding of 
these newly-identified pulsars. First, one gets a 
measurement of the period and period derivative 
of the source. Having these two numbers allows us 
to derive estimates of several key parameters in- 
cluding the characteristic age, the inferred dipole 
magnetic field strength, and the spindown energy 
loss rate. These parameters are fundamental to 
understanding the astrophysics of the system. For 
example, the characteristic age is useful (though 
certainly not definitive) in the context of argu- 
ments for or against associations with supernovae 
and pulsar wind nebulae (PWNe). 

The next critical parameter in the timing model 
is the pulsar position. Estimating the source posi- 
tion from the reconstructed photon arrival direc- 
tions can yield localizations that are good to a few 
arcmin, but to do better than this requires timing. 
For young or middle-aged pulsars, the LAT can 
measure pulse arrival times with accuracies of or- 
der a millisecond^ which can be fit to determine 
positions to arcsecond accuracy. Accurate posi- 
tions then allow deep counterpart searches in the 
X-ray, optical, and radio bands, and remove the 
effects of position error on the remaining timing 
parameters, most notably the spin-down rate. 

Once the basic spin and position parameters 
are well determined, timing allows us to inves- 
tigate the rotational irregularities that are com- 
mon in young pulsars. The primary phenomena 
are timing noise and glitches. Glitches are sud- 
den increases in pulse frequency with a magnitude 
in the range /S.v/v ~ 10~ 10 — 10~ 5 , which pro- 
vide valuable information about the superfluid in- 
terior of neutron stars ( Andersson et al.|2003 Link 



^The accuracy of a pulse time of arrival measurement is 



determined by the photon statistics and the sharpness of 
the features in the pulse profile. It is always considerably 
larger than the ~ 1 /is accuracy on individual photon event 
times recorded by the LAT. 
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et al.|[l999 for example). Timing noise is unmod- 
eled low-frequency (often quasi-periodic) noise ob- 
served in the residuals of many pulsars after all 
the deterministic spin-down effects have been re- 
moved. The magnitude of the timing noise has 
been shown to correlate with frequency derivative 



(i.e. torque) ( Cordes & Downs 


1985 


Arzoumanian 


et al.|1994a |Hobbs et al.|2010 


1, but its nature re- 



mains poorly understood. 

Using the timing positions for these pulsars, we 
have also undertaken deep radio observations of 
the 7-ray-selected pulsars to search for radio pul- 
sations. These searches have resulted in three dis- 
coveries of radio pulsations, which have been pub- 



lished elsewhere (Camilo et al. 2009 Abdo et al 



2010c I. Here, we compile the upper limits from 
our observations, and from the literature, for the 
remaining pulsars. These are important inputs to 
population statistics and modeling of these appar- 



ently radio quiet pulsars (Yadigaroglu & Romani 
1995 Story et al.|2 007 for example) , as well as for 
guiding future deeper searches. 

2. Methods 

2.1. Data Selection 

For the current analysis we use LAT data from 
2008 August 4 through at least 2010 February 4, 
the first 18 months of LAT survey operations. We 
select LAT events from the most restrictive "dif- 
fuse" class of the "Pass 6" event reconstructions 
( |Atwood et al.|2009[ ) with a zenith angle of < 105° 
to reduce contamination from atmospheric sec- 
ondary 7 rays from near the Earth's limb. For each 
pulsar, we find an optimal radius and low energy 
cut to maximize the pulse detection significance. 
The radius cuts ranged from 0.5-1.6°, while the 
low energy cuts ranged from 50-900 MeV. Only 
photons that pass these cuts are included in the 
timing analysis. 

The number of photons surviving these cuts 
ranged from 1,174 (PSR J0633+0632) to 14,875 
(PSR J1836+5925) in our 18 months of observ- 
ing, a span in which the pulsars completed several 
hundred million rotations. This emphasizes the 
unique nature of timing pulsars using extremely 
sparse 7-ray data. Typically only of order 100 
photons go into each TOA determination. In addi- 
tion, unlike with radio pulsar timing, the integra- 
tion time per TOA is equal to the spacing between 



TOAs, requiring the model to maintain phase ac- 
curacy over a much longer time than is required for 
radio pulsar timing where the integration time for 
a TOA is only minutes or hours. We constructed 
initial models using the prepf old tool from the 
PRESTO pulsar analysis software packag^J This 
tool performs epoch folding searches over narrow 
ranges of frequency and frequency first and sec- 
ond derivatives around the v and v values found 
from the blind search to maximize the signal-to- 
noise ratio. Combined with searching over a grid 
of possible pulsar positions, we are able to arrive 
at an initial model that maintains coherence well 
enough for TOAs to be determined and the pulsar 
timing to proceed as described below. 

2.2. Geocentering 

Pulsar timing software generally expects pulse 
TOAs to be measured at an observatory that is at 
a fixed geographic location on the Earth. Obser- 
vations from a spacecraft in orbit about the Earth 
obviously do not satisfy this condition and this 
must be accounted for before computing a pulse 
arrival time. One could go directly to a time scale 
(such as TDB) at the solar system barycenter, but 
this requires a precise knowledge of the pulsar lo- 
cation before the correction can be done and re- 
moves the possibility of fitting for the pulsar po- 
sition as part of the timing model. Instead, in 
order to remove the effects of the spacecraft mo- 
tion on the photon arrival times while maintaining 
the ability to fit for astrometric parameters in the 
timing model, we correct the measured times to a 
fictitious observatory located at the Earth's geo- 
center. 

LAT photon times are recorded in Mission 
Elapsed Time (MET), which is referenced to Ter- 
restrial Time (TT) via the MJDREF keyword 
111 the FITS file headed Time is maintained on- 
board the spacecraft to an accuracy of better than 
1 [i% using a GPS receiver dSmith et al.1[2008| . 

The geocentric time is the satellite time cor- 
rected for geometric light travel time to the geo- 
center. It does not include relativistic terms in 



from Ihttp: //www. cv.nrao . edu/~sransom/ 



2 Available 
presto/ 

3 See OGIP Memo OGIP/93-003 http://heasarc.gsfc 



nasa . gov/docs/heas arc/of wg/docs/rates/ogip_93_003/ 
ogip_93_ 003 ■ html| 
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the correction. The geocentric photon time t geo is 
defined as 



t 



gco 



to 



■bs 



n 



(1) 



where r sat is the vector pointing from the geocen- 
ter to the spacecraft, n psr is a unit vector pointing 
in the direction of the pulsar (here assumed to be 
at an infinite distance), and c is the speed of light. 

This correction is applied using the Fermi sci- 
ence tooQ gtbary with the tcorrect=geo op- 
tion. After this correction, the time system for 
the events is still TT, but all times are referenced 
to the geocenter. This correction has a maximum 
amplitude of 23.2 ms. Therefore, an error in the 
assumed pulsar direction as large as 1° causes a 
maximum error in the corrected time of only 0.4 
ms. 

2.3. TOA Determination 

A TOA is determined from the photon times in 
a segment of data by first assigning pulse phases 
to each photon based on an initial model, then 
measuring the phase offset (A) required to align a 
standard template profile with the measured pulse 
profile (see Figure]!]). This offset is then converted 
to a time using the pulse period (P) and added 
to the observation start time, To, to become the 
measured TOA. 



TOA = T + A x P 



(2) 



This TOA is the time when the fiducial point on 
the pulse profile arrived at the observatory, for a 
representative pulse during the observation inter- 
val. In the case of the geocentered LAT events 
the TOA is for a fictitious observatory at the geo- 



center (observatory code coe in Tempo2; Hobbs 
et al. ( |2006 >). The measurement can be made us- 
ing binned pulse profiles (as in Figure[T]) or directly 
from the unbinned photon phases, as described be- 
low. 

For this work, we divide the full 18-month ob- 
servation interval into segments of equal duration 
and determine a pulse time of arrival from each 
segment. The length of each segment is a balance 
between signal-to-noise ratio and time resolution. 
Longer integrations result in better signal-to-noise 



http : //f ermi . gsfc .nasa. gov/ssc/dat a/analysis/ 



documentation/ 



ratio and smaller statistical measurement errors 
on each TOA. On the other hand, shorter inte- 
grations provide finer time resolution that better 
samples the annual sinusoidal signal caused by the 
Earth's motion around the Sun and the timing 
noise in some very noisy young pulsars. There- 
fore we try to achieve at least 1 TOA per month. 
Fainter pulsars that require a substantial fraction 
of a year, or longer, per TOA measurement will 
be difficult to time with the LAT. 

For most radio pulsar timing, the TOAs are 
determined from binned data. The start time 
of the observation is precisely known from the 
observatory clock. During the observation, data 
are folded using predicted phases for the pulsar 
based on a provisional ephemeris (e.g. using the 
-polyco option to Tempo2), and a binned pro- 
file for that observation is computed. The arrival 
time is computed by cross-correlating the observed 
profile with a high signal-to-noise template pro- 
file with the same binning. The accuracy of this 
measurement is improved if the cross-correlation 
is implemented as a fit to a linear phase gradi- 
ent in the Fourier domain (an application of the 
Fourier shift theorem), rather than as a simple 
time-domain cross correlation (Taylor 1992). Fi- 



nally the TOA is determined as the observation 
start time plus the measured phase offset (con- 
verted back into time units). 

The binned TOA determination method can 
also be applied to photon data, such as that from 
the LAT, by computing the predicted phase for 
each photon and building a binned pulse profile 
from the events. However, since we must make 
TOA measurements based on a small number of 
detected photons (often < 100 photons go into 
each TOA), we can improve the TOA determi- 
nations by using an unbinned likelihood analysis 
to compute the TOA directly from the set of pho- 
ton phases. This has been discussed before (Liv- 



ingstone et al. 2009), but we have developed and 



generalized the technique and describe it in detail 
here. 

We use an unbinned maximum likelihood 
method to estimate both the light curve tem- 
plate and the TOAs with associated errors. In 
the likelihood formulation, the template is inter- 
preted as a periodic probability density function 
to observe a photon at a given phase, /(</>; A, A), 
with A some set of parameters describing the light 
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curve morphology and A accounting for the phase 
shift between the template and the given data set. 
The TO A is determined by A. The template is 
normalized such that J d<f> f((f>; A, A) = 1. 

We first start with a description of the tem- 
plate, /(</>; A, A), which must be a continuous func- 
tion that can be evaluated at any value of 0. In 
many cases, the statistics arc sufficiently limited 
that the pulse profile can be described as a sum 
of a constant background component and a small 
number of gaussian peaks. That is, 



the coefficients for the fcth harmonic are 



f{<f>) = (!- ^ZPi) 9{<t>,x % ,a % ), (3) 

i=l i=l 

with N p the number of gaussian peaks, pi the frac- 
tion of the total emission belonging to each peak, 
and g(x, a) a gaussian with mean x and standard 
deviation a. The domain of the gaussian functions 
is assumed to be wrapped to [0, 1). Here, A can be 
associated with x\, the location of the first peak, 
while the remaining parameters are subsumed in 
A. 

With increasing statistics, the complexity of 
GeV light curves is no longer well-represented 
by a simple sum of components. Bridge emis- 
sion and peak asymmetry appears, and in gen- 
eral, no simple functional form is sufficient to de- 
scribe the profile (e.g. the Vela pulsar (Abdo 



et al. 2010el). In this case, we prefer kernel den- 
sity estimation (KDE) methods (de Jager et al. 
1986). These methods result in a faithful, non- 



parametric representation of the light curve. How- 
ever, even for bright pulsars, the available statis- 
tics are such that KDE methods produce a tem- 
plate with broadened peaks and "noisy" valleys, 
neither of which is desirable for the calculation of 
a TOA. A good estimator for the template should 
provide a smooth template (ignore fluctuations) 
while simultaneously preserving the structure and 
sharpness of the peaks, which is important because 
the template sharpness is a factor in the accuracy 
of the TOA measurements. We outline two ap- 
proaches we have found effective below. 

The first forms the basis for the H-test statis- 
tic often used to assess pulsation significance in 



the absence of a template (de Jager et al. 19891 



Coefficients of a Fourier expansion are estimated 
directly from the unbinned phases. For n photons, 



Ecos(fc<fo), A = - } sin(k(f>i 



(4) 



i=l 



i=l 



and the light curve is given by 



f{4>) = l + 2j2^k cos(£^) + fa wx(hf>). (5) 

k=l 

The only free parameter is the overall phase of the 
light curve; variation can be implemented with the 
Fourier shift theorem or simply by adding a con- 
stant phase to the data. The number of harmonics 
retained should offer an optimum balance between 
peak "sharpness" and noise in the remainder of the 
profile. We call this the 'empirical Fourier' (EF) 
method. 

The second method is a gaussian KDE with a 
phase- dependent bandwidth, the idea being to use 
smaller bandwidth for the peaks while smooth- 
ing the valleys with a broader kernel. Here, 

/(<£) = S"=iff( ( ^ < ^> cr i): witn 9 again the stan- 
dard gaussian. The bandwidth is determined by 
c>i = (/max - /min)//(0i)\Ai- Lest the reader 
worry about this circular definition, in practice 
we begin with a phase-independent bandwidth 
a = y/n and iterate. As with the Fourier expan- 
sion, the only free parameter is the overall tem- 
plate offset. 

In summary, for our template, we choose one 
of the three above methods (Gaussian, EF, KDE) 
that produces the best results, as evidenced by the 
smallest RMS residuals. A comparison of a pulse 
profile fitted with the three different templates is 
shown in Figure [2| The template choice is docu- 
mented for each pulsar. 

With the template defined, the next step is to fit 
for the TOA from each segment of data, always us- 
ing the chosen template to define the pulse profile 
and fiducial point. For the fitting, we start with 
an approximate timing solution (say from a v-v 
search) and fold the photon arrival times to obtain 
a set of phases {<f>i, <f>2, 4> n }- The probability to 
observe these data, given the light curve model, 
is formally inverted to form the log likelihood for 
the parameters, J27=i 1°§ /(^> ^> 4>i)- The param- 
eters are varied to maximize the log likelihood. 
In the case of a multi-gaussian template, the full 
dataset is used to determine A, while in determin- 
ing TOAs, only A is fit. The likelihood surface 
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generally has a gaussian shape near the best-fit 
value for A, and we estimate the error on A by 
measuring and inverting the curvature of the log 
likelihood function at the best-fit value. Thus, to 
determine a TOA, the above fit is carried out for 
each subset to estimate for A and its error a&. 

We mention here an additional challenge brought 
on by morphology of GeV light curves. An appre- 
ciable fraction of pulsars observed so far present 
light curves with two peaks of similar height with a 
separation close to 0.5 periods. These light curves 
are approximately invariant under a half-period 
translation, and statistical fluctuations may then 
lead to a likelihood maximum associated with the 
"wrong" peak. When this happens, a blind search 
for the maximum likelihood will result in a TOA 
off by 0.5/is seconds which must be excluded from 
the timing solution fit. To avoid loss of data, 
rather than employing a blind search, we "track" 
the solution. That is, provided the trial solution is 
sufficiently good (and this is always the case with 
iteration), the drift of the actual arrival time from 
the predicted arrival time is much less than half 
of a period. We then simply restrict the search for 
the likelihood maximum to within a range that 
excludes the "wrong" peak. 

2.4. Fitting Timing Models 

The measured sets of TOAs are then fitted to 
a timing model using the pulsar timing software 
Tempo2 flHobbs et al."p006| [Edwards et al.|2006 |. 
There arc many parameters that can be used in the 
timing models. For all pulsars we fit for pulse fre- 
quency (v) and frequency first derivative (z>) and 
frequency second derivative (i>). We fit for v as a 
measure of the timing noise present in each pulsar. 
In most cases, a significant v is not detected and 
we report a 2 a upper limit on the magnitude \v\. 
In the cases where v is measured, we attribute this 
solely to timing noise as the v expected from any 
reasonable braking index would be immeasurable 
over our 18 month data span. If this is still in- 
sufficient to whiten the residuals, we add a third 
frequency derivative, or harmonically related si- 



nusoids (WAVE parameters in Tempo2, see Hobbs 
et al. ( |2004 )) to the fit until a satisfactory model is 
achieved. Note that because of large covariances 
between the parameters, one should avoid fitting 
the position and WAVE parameters at the same 
time. 



Finally, in three cases (PSRs J0007+7303, 
J1124-5916, J1813-2332), a glitch was observed 
and several glitch parameters were added to the 
fit, as described in <|3] 

In our models, the absolute phase 0.0 is arbi- 
trary. In the case of pulsars with both radio and 
7-ray emission the convention is usually to assign 
phase as the peak of the radio pulse (as a proxy 
for the more physically meaningful point of closest 
approach of the magnetic axis to the line of sight). 
However, since we don't observe radio pulsations 
from most of these pulsars we have not attempted 
to define a particular phase 0. However, we do 
report the parameter TZRMJD for our models, 
which is the reference for phase 0.0. Phase 0.0 is 
the pulse phase at the time TZRMJD at the geo- 
center at infinite frequency. 

As a final note, we want to emphasize that dif- 
ferent timing models are appropriate for different 
purposes. One of the primary goals of this work is 
to use the capability of pulsar timing with the LAT 
to make accurate localizations of 7-ray-selected 
pulsars, thus enabling multiwavelength studies of 
potential counterparts. A secondary goal is char- 
acterizing the timing noise in this set of pulsars. 
For other purposes, different models are appro- 
priate. In particular, for many studies it would 
be preferable to freeze the position using an ac- 
curately known counterpart position (from Chan- 
dra X-ray observations, for example) to reduce the 
number of free parameters in the model. 

The Tempo2 timing models described here will 
all be made available electronically at the Fermi 
Science Support Center (FSSC) web sit<^] 

2.5. A Tempo2 Plugin For Assigning Pho- 
ton Pulse Phases 

One important use of the timing models pre- 
sented here is to be able to assign an accurate 
pulse phase to each photon in a LAT observation 
of a particular pulsar. This is needed for stud- 
ies of the 7-ray light curve, phase-resolved spec- 
troscopy, or "gating" the data on the off-pulse re- 
gion to blank out a pulsar to enable studies of faint 
sources nearby (e.g. Cyg X-3 ( Abdo et al.|2009c )). 
The standard Fermi Science Tool gtpphase tool 
was developed for this application. However it 



'http: / /fermi. gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
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suffers from the limitation that it cannot repre- 
sent the full complexity of pulsar timing mod- 
els that include frequency derivatives above i>, 
glitches, parallax, proper motion, or WAVE param- 
eters. Tempo2, on the other hand, allows all of 
these as well as several additional orbital models 
for pulsars in binary systems. 

For these reasons, we have implemented a 
graphical plugin for calculating pulsar phases for 
Fermi-L AT data with Tempo2, called f ermi.plug.C. 
This plugin takes LAT event ("FT1") files with 
the photon arrival dates, spacecraft ("FT2") files 
with the satellite position as a function of time, 
and Tempo2 timing solutions and writes photon 
pulse phases in the FT1 event file. It uses the 
same spacecraft position interpolation algorithm 
as implemented in the Fermi science tool gtbarjj^] 
and derives barycentric photon times with anal- 
ogous methods. The barycentric times are then 
treated as TOAs to find the pulsar phases relative 
to the absolute phase reference given by the TZR- 
MJD parameter in the input ephemeris. The plu- 
gin thus allows Fermi-LAT data analysis with an 
ephemeris built from radio, X-ray or 7-ray TOAs, 
with virtually unlimited complexity in the timing 
model. The plugin has been shown to reproduce 
the results from the Science Tools when working 
in tempol emulation mode, but using this mode 
is not required. It is available in the Tempo2 
sourceforge distribution^] and from the Fermi Sci- 
ence Support Centci]^] This plugin is suitable for 
use with any of the timing models presented here. 

3. Results 

In this section, we present details of the timing 
models for each of the 17 pulsars listed in Table 
[TJ The models arc determined from the data set 
as described above. The statistical errors on the 
parameters are the (single parameter l-er) uncer- 
tainties reported by Tempo2 from the fits. For the 
pulsars with no v required in the model, we esti- 
mate the statistical error in the position fit from a 
fit with v free, because this results in a more con- 
servative error estimate that better accounts for 
the correlations between the astrometric and spin 



6 See http : / /f ermi . gsf c . nasa.gov/ssc/library/ support/ 



psr_tools_anatomy/ 

7 http:/ /tempo2. sourceforge. net/ 



^http : //f ermi ■ gsf c ■ nasa. gov/ssc/data/ analysis/user/ 



parameters. 

The errors on the position in the tables and 
shown in the figures are statistical only (though in 
the figures they are 95% confidence, rather than 
1-cr since that is the standard practice for LAT er- 
ror ellipses), and thus underestimate the true error 
on the position determinations. This is predomi- 
nately because our span include only 1.5 periods 
of the annual sinusoid induced by an error in the 
position. Timing instabilities present on a similar 
timescale can thus perturb the fitted position. In 
addition, there arc strong covarianccs between the 
astrometric and spin parameters that mean that 
the parameters are not actually determined to as 
high a precision as indicated by the 1-parameter 
statistical errors. Because knowing the true posi- 
tional uncertainties is very important for counter- 
part searches at other wavelengths, we have tried 
to quantify the magnitude of this effect by Monte 
Carlo simulation. We make the assumption that 
over the short span of data we have, the measured 
values of v and v are dominated by timing noise, 
not the secular spindown of the pulsar. To esti- 
mate the magnitude of the systematic error, we 
generate many simulated sets of TOAs, each one 
using the measured timing parameters for the pul- 
sar, with the exception of v and v . For those 
two parameters, we replace them with normally- 
distributed random values with mean zero and 
standard deviation equal to the measured value, 
or the upper limit in the case where v is not signifi- 
cantly detected. For v and v we use random values 
distributed around the measured value with the 
measured uncertainty. Each trial set of TOAs is 
then fit with a 1-year sine wave plus a polynomial 
up to order v (see Figure [3]) and the magnitude of 
the sine wave is converted to a position offset (see 
Appendix [A]) . We then compute the position un- 
certainty by finding the position offset that 68% 
of the trials are lower than. It is important to 
note that our simulations include random, uncor- 
related, measurement errors as appropriate for the 
particular pulsar, so these position error estimates 
are of the total uncertainty, including both system- 
atic and statistical components. Also, the fidelity 
of the estimates depends on how well our random 
polynomial model of the timing noise describes the 
actual situation, which is not well understood and 
may vary from pulsar to pulsar. Therefore, these 
estimates should be considered indicative of the 
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magnitude of the position error, but not precise 
bounds on the systematic errors. Finally, in this 
analysis, we just consider the total position off- 
set, so it yields an intermediate value in the cases 
where the position error region from the timing is 
highly elliptical. 

For all timing models we use the JPL DE405 
planetary ephemeris (IStandish 1998). All reported 



frequencies and epochs are referenced to the TDB 

as has been the 



time system (Seidelmann 



1992) 



standard for pulsar wor The clock correction 
procedure is TT(TAI) and all fits are made with 
weighting by the TOA error estimates enabled 
(MODE 1). The reference time (TZRMJD) is for 
the geocenter at infinite frequency. The valid- 
ity range for each model is included in the ta- 
bles. Care should be taken when attempting to 
use these models outside of that range. In par- 
ticular, models that include significant v or higher 
order derivatives, or WAVE parameters, will extrap- 
olate very poorly outside the fit range, since timing 
noise is a stochastic process and neither of those 
parameterizations reflect a physical model of the 
process. 

In the following subsections, we discuss some 
of the main results from the timing of the individ- 
ual pulsars. For each pulsar, we present the timing 
model in a table, the post-fit timing residuals, a 2- 
D phaseogram, and a pulse profile. In all cases, the 
optimized data selections used for the pulsar tim- 
ing are also used to construct the 2-D phaseogram 
and pulse profile figures. The phaseograms are raw 
photon counts and are not exposure-corrected, so 
the apparent variations in brightness that can be 
seen are from exposure variations resulting from 
the ~ 55-day precession period of the spacecraft 
orbit, the change in rocking angle during the mis- 
sion, spacecraft reboots, or automatic repoints in 
response to 7-ray bursts. The fluxes of 7-ray pul- 
sars is expected to be constant on time scales of 
days to months. 



In the discovery paper (Abdo et al. 2009a I, 
five of the pulsars were given names of the 
form JHHMM-I-DD because the declinations were 
not known with sufficient precision to justify a 
name of the form JHHMM+DDMM. In all five 



cases, we now know the position well enough to 
confidently add the additional precision to the 
names, as shown in Table [T] Also, in five cases 
(PSRs J1418-6058, J1741-2054, J1809-2332, 
J1813-1246, and J1958+2846) the current best- 
fit timing position would result in a different last 
two digits of the declination than given in the dis- 
covery paper, although in several cases we know 
the name to be correct based on the X-ray coun- 
terpart position. In all cases, we follow the IAU 
preference for not changing a source name once it 
is given and we continue to use the original names, 
except where we have only added precision, as de- 
scribed above. See the sections on each individual 
source below for a discussion of the confidence in 
the previously proposed counterpart associations. 

3.1. PSR J0007+7303 

The timing model parameters for this pulsar 
are displayed in Table [4] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [5] |6j and [7j respectively. 

This was the first pulsar discovered in a blind 



search of 7-ray data (Abdo et al. 2008) and is 



believed to be the pulsar powering the compact 
PWN RX J0007.0+7303 near the center of the 
shell-type supernova remnant CTAl. As seen in 
Figure [5j our timing position provides indepen- 
dent confirmation of that conclusion. In addition, 
we have detected a glitch in this pulsar on 2009 
May 1 with a magnitude Avjv = 5.53(1) x 10~ 7 , 
a typical glitch magnitude for a pulsar of this age. 
When we fit for position in the timing model, the 
glitch can be fully accounted for by a simple Av 
at the time of the glitch. However, when we hold 
the position fixed at the Chandra position of the 



9 Note that the default time system for Tempo2 is TCB, but 
we override this default and use TDB as the time units for 
our models. 



point source (00:07:01.56, 73:03:08.3; see |Halpern 
eTaLlpoMl )) we find that an additional parameter 
is required. This can be modeled as a change in 
the frequency first derivative at the glitch of Avjv 
of 0.0010(2). It is important to note that v and 
the glitch Av are highly covariant and additional 
data will likely be required to determine whether 
timing noise or a frequency derivative change at 
the glitch are the correct model for the observed 
behavior. The properties of this source and the 
glitch will be discussed in more detail in a future 
paper (Abdo et al. 2011, in prep). 
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Table 1 

Pulsars Timed with the Fermi LAT 



JName 


Frev. IName 


Fenod 








(ms) 


( 1 n34 . ^ — 1 ^ 

(10 erg s J 


J0007+7303 




315.9 


45.2 


J0357+3205 


J0357+32 


444.1 


0.6 


J0633+0632 




297.4 


11.9 


J1124-5916 




135.5 


1195.0 






1 1 n fi 

1 ±u .u 


4Q4 8 

lyi . O 


J1459-6053 


J1459-60 


103.2 


90.9 


J1732-3131 


J1732-31 


196.5 


14.5 


J1741-2054 




413.7 


0.9 


J1809-2332 




146.8 


42.9 


J1813-1246 




48.1 


624.1 


J1826-1256 




110.2 


358.0 


J1836+5925 




173.3 


1.1 


J1907+0602 


J1907+06 


106.6 


282.7 


J1958+2846 




290.0 


34.2 


J2021+4026 




265.3 


11.6 


J2032+4127 




143.2 


27.3 


J2238+5903 


J2238+59 


162.7 


88.9 



3.2. PSR J0357+3205 




Fig. 1. — Example of a TOA measurement. The 
blue histogram is a binned pulse profile generated 
from the observed photons (two cycles are shown 
for clarity) . The red curve is a two Gaussian tem- 
plate profile, where the point at phase 0.0 (or 
equivalently 1.0) defines the fiducial point. The 
black arrow represents the measured phase offset 
(A) required to align the profile with the template. 



The timing model parameters for this pulsar 
are displayed in Table [5] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [8j |9j and [lOj respectively. 

PSR J0357+3205 is the slowest spin period (444 
ms), and lowest E (5.8 x 10 33 erg s _1 ) pulsar in 
our sample. In the discovery paper (Abdo et al 



2009a), it was flagged as having a potentially large 



systematic error in the i> and the parameters de- 
rived from it, because of the uncertain position. 
The long period, low count rate, and relatively 
broad pulse profile still limit the timing precision 
to an RMS of 5.3 ms, but nevertheless the fre- 
quency derivative is now determined to an accu- 
racy of ~ 0.2 percent. 

For this low E, the distance is constrained to 
be < 870 pc, assuming the flux correction factor 
/n = 1 ( |Watters et~aLl|2009[) and using t he LAT 
7-ray flux (Gioo) from Abdo et al. ( 2010d[ ) to keep 



the 7-ray efficiency < 1. As seen in Figure |5J no 
X-ray counterpart is apparent in a Swift image 
of the region, which is not surprising in such a 
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1.0 

Pulse Phase 



Fig. 2.— The pulse profile of PSR J0007+7303 fit- 
ted with a 2 Gaussian, Kernel Density Estimator 
(KDE), and Empirical Fourier (EF) template with 
16 harmonics. The black histogram shows the 
measured pulse profile with 32 bins, but the tem- 
plates are fitted to the unbinncd photon phases, 
as described in the text 



shallow exposure. However, as the pulsar is at 
such a small distance, this is a promising target 
for deeper XMM-Newton or Chandra follow up. 
Using the Tempo2 simulation capability, we pre- 
dict a 1-ct uncertainty on the timing position of 2" 
after 5 years of observation. 

3.3. PSR J0633+0632 

The timing model parameters for this pulsar 
are displayed in Table [6] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures (TTJ [l2j and [l3j respectively. 

This pulsar is also rather slow (297 ms) and 
faint (only 815 photons detected per year), but 
the timing is still quite good (RMS = 1.4 ms), 
as a result of the very narrow pulses. The tim- 
ing localization is close to the X-ray point source 
Swift J063343. 8+063223 which was proposed as 



the counterpart by Abdo et al. (2009a 



To further study the X-ray counterpart, we ob- 
tained a 20 ks Chandra ACIS-S image of the region 
on 2009 December 11 (ObsID 11123). The X-ray 
point source counterpart to the pulsar is clearly 
visible in Figure 11 We measure a position of 
06:33:44.142, +06:32:30.40, which is 4.6" from the 
best fit timing position. To fit the spectrum of this 
source, we analyzed the data using CIAO version 
4.3 with the latest calibrations (CALDB 4.1.1) 
applying the standard particle background sub- 
traction and exposure correction. We extracted 
326 photons from a 3.5 pixel extraction region 
around the source location (for a count rate of 
1.63 x 10~ 2 cts s . We see no evidence for a 
compact (arcsecond-scale) PWN in the immedi- 
ate vicinity of the point source. To fit the spec- 
trum, we found that an absorbed blackbody + 
powerlaw model is required. We obtain the follow- 
ing parameters from our fits, with 90% confidence 
error estimates: n# = 0.15to'io x 10 22 cm~ 2 , 

This model 



kT = 0.1lto;°2 keV, r = 1.5 ± 0.6 



yields a 0.5-8 keV flux estimate of 9.2+^ x 10 



-14 



If we instead fit an absorbed neu- 



erg cm t 

tron star atmosphere (nsa; Zavlin et al. (1996)) 



plus powerlaw model, we find a somewhat higher 
tih of 0.24^Q 2 2 x 10 22 cm -2 , a lower temperature 
of kT = 0.048^;^ keV, and a similar photon 
index T = l-39±g;l- 

To look for larger-scale extended emission, we 
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smoothed the Chandra image with a gaussian ker- 
nel of 1.5" width (see Figure [4]) and find a faint 
X-ray nebula extending about an arcminute south 
of the pulsar. In the region of the PWN (as shown 
in the Figure [4]), we find an excess of 738 counts 
on a background of about 1600 counts and have fit 
the integrated spectrum with an absorbed pow- 
erlaw model. With all parameters free, we find 
n H = 0.1±g:f x 10 22 cm" 2 , T = 0.9t£f for a 
flux in the 0.5-8 keV band of 2.2 ± 0.5 x 10" 13 
erg cm -2 s _1 , where the error regions are at the 
90% confidence level. If instead, we freeze uh at 
0.154 x 10 22 cm -2 , as found in the blackbody + 
powerlaw spectral fits of the point source, we find 
a 90% confidence range for the photon spectral 
index T of 0.74-1.29. 

3.4. PSR J1124 5916 

The timing model parameters for this pulsar 
are displayed in Table [7] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [Mj [l5j and [THJ respectively. 

This pulsar with very small characteristic age 
(r c = P/2P = 2900 yr) is associated with the 
supernova remnant G292.0+1.8 and is the only 
one in this sample that was previously known as 

has a 



al. 



2002). 

T 



It 



a radio pulsar (Camilo et 

very high E of 1.2 xl0 aY erg s _i and exhibits 
a great deal of timing noise. It is also is very 
faint, with a 1.4 GHz flux density of only 0.08 



mJy (Camilo et al. 2002) and far enough south 



that it can only be timed with the Parkes Tele- 
scope, where it requires several hours of integra- 
tion to even get a detection. Therefore, it has 
not been regularly timed with radio observations 
since its discovery. Because there was no contem- 
poraneous radio ephemeris available, the LAT pul- 
sations from this source were discovered using a 
limited blind search around the known spin pa- 
rameters. With our LAT timing, we are able to 
obtain TOA uncertainties of 1.5-3.2 ms every two 
weeks and obtain a phase-connected timing model. 
The timing position is 1.8" from the Chandra 
source CXOU J112439. 1-591620 flCamilo et al 



2002). The pulsar exhibited a glitch of magni- 
tude Avjv = 1.6 x 10~ 8 around MJD 55191. A 
small Av/v of —0.00472(3) was also observed at 
the glitch. 

The very large measured v results in a Monte 



°'^|600 54700 54800 54900 55000 55100 55200 55300 
Time (MJD) 



Fig. 3. — Example of a Monte Carlo trial to esti- 
mate the systematic error on the timing position. 
The blue points with error bars are the simulated 
TOAs, which are fitted to the sum of a 1-year sinu- 
soid (solid red line) and a third order polynomial 
(dashed green line), as described in the text. 




51.00s 48.00s 45.00s 42.00s 

RA (J2000) 



6h33m39.00s 



Fig. 4.— Chandra 0.5-7 keV ACIS-S image of 
PSR J0633+0632, smoothed with a gaussian of 3 
pixels (at a scale of 0.5" per pixel) to highlight 
the extended PWN emission. The extraction re- 
gion used for the PWN spectral analysis is shown, 
where the two point sources are excluded from the 
region. The pulsar is the northernmost of the two 
point sources in the region. 
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Carlo position error estimate of about 1', but 
given the good agreement between the timing po- 
sition and the Chandra position, this must be a 
large overestimate, perhaps because the assump- 
tions inherent in the Monte Carlo estimate are vi- 
olated. The measured v results in a braking index 
n = = —3.78, which is of comparable magni- 
tude, but opposite in sign to the n = 3 expected for 
vacuum dipole braking ( Lorimer fc Kramer|2005 1 . 
This measurement, along with the substantial red 
noise still present in the timing residuals (see Fig- 
ure 15 1 all suggest that the spindown of this pulsar 



is rather noisy. 

To measure the radio to 7-ray phase alignment 
of this pulsar, we made a 5-hour observation with 
the Parkes Radio Telescope at a frequency of 1.4 
GHz. Since this pulsar is not timed routinely in 
the radio, we required a new, contemporaneous, 
observation because the extreme timing noise in 
this system prevents the timing model from be- 
ing extrapolated forward or backwards in time. 
This radio light curve was presented previously in 
the First Fermi LAT Catalog of Gamma-ray Pul- 



sars (Abdo et al. 2010d), but the absolute phase 



alignment presented there was incorrect. The ver- 
sion presented here correctly accounts for the de- 
lay from interstellar dispersion using DM = 330 
pc cm -3 . The new value for the lag from the ra- 
dio peak to the first 7-ray peak (6) is 0.128(3). 
This correction was also made to the catalog in an 
erratum ( |Abdo et aL]|2011 ). 



3.5. PSR J1418-6058 

The timing model parameters for this pulsar 
are displayed in Table [8] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [l7j [l8j and [l9j respectively. 

In the discovery paper this pulsar was proposed 
to be associated with the PWN G313. 3+0.1 ("The 



Rabbit"). Ng et al. (20051 find two point sources 



(Rl and R2) in Chandra and XMM observations 



of the region. Recently, Roberts ( 2009 ) reported a 



weak detection of X-ray pulsations at the period of 
PSR J1418-6058 in XMM data, so it is believed 
that Rl is the correct counterpart. However, the 
timing position is 14" from the position of Rl (RA. 
= 14:18:42.7, Decl. = -60:58:03), which is signif- 
icantly larger than the 2" statistical error on the 
position. It is important to note that this pulsar is 



very noisy and the timing model, which includes 
terms up to the frequency third derivative, clearly 
does not fully describe the data (see Figures 18 
and 19 1. This causes the statistical error to be 



underestimated and there is a significant system- 
atic error on the position as well. For example, 
adding a fourth frequency derivative to the model 
causes the position to shift by 10". As seen in Fig- 
ure |17| there are three X-ray point sources near the 
nominal timing position, but none are coincident 
with the timing position. The X-ray point source 
R2 is not apparent in this image, so it is likely 
variable, and probably not the counterpart to the 
pulsar. A Monte Carlo estimation of the system- 
atic error (see Sj3| on the position induced by the 
timing noise seen in this pulsar is 7" for the poly- 
nomial model for timing noise and 40" using the 
red noise model with an RMS of 62 ms (as com- 
puted from the measured v and v for the pulsar) . 
Therefore, we can't exclude Rl as being the coun- 
terpart based on the positional disagreement. In 
addition, the faint X-ray source just north of the 
pulsar, which we call CXOU J141843.3-605734, is 
equally consistent with the timing and so further 
data, or a confirmation of the X-ray pulsations 
from Rl, will be required to confirm the associa- 
tion with either source. 

3.6. PSR J1459 6053 

The timing model parameters for this pulsar 
are displayed in Table [9] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [20j [21] and [22] respectively. 

In the discovery paper, no counterpart was 
known for this pulsar. With the high precision 
timing now available, we see that the Swift image 
shows an apparent faint point source near (9.8" 
offset from) the timing position (see Figure [20| . 
We call this source Swift J145931. 3-605319, but 
its properties are not well constrained because of 
the faintness in the 6 ks Swift image. A deeper 
X-ray image is required to confirm this source. 

3.7. PSR J1732 3131 

The timing model parameters for this pulsar 
are displayed in Table [10] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
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in Figures 23] 24 and 25 respectively. 

This source shows minimal timing noise, with 
only an upper limit of 2 x 10~ 23 s~ 3 on The 
timing error ellipse is significantly elongated be- 
cause of the low ecliptic latitude of the source 
(/? = —8.2°), causing the declination to be more 
poorly constrained than the R.A. Earlier Swift 
imaging showed no significantly-detected X-ray 
source at the pulsar location, so we pursued a 
deeper observations with Chandra. Our 20 ks 
Chandra ACIS-S image (ObsID 11125) reveals an 
X-ray point source consistent with the timing po- 



sition (see Figure 23). We measure the position as 
17:32:33.551, -31:31:23.92, which is 0.9" from the 
timing position, well within the 95% confidence 
region. 

We performed a spectral analysis of the source 
based on 79 photons from the point source (with 
< 1 count from the background). Since the source 
is still detected even with a low energy cut of 
3.5 keV, it is clear that a non-thermal compo- 
nent is required. However, with the small number 
of counts, the power law photon index cannot be 
constrained, so we freeze it at V — 1.5 in the fits. 
The spectral parameters from our fits to an ab- 
sorbed blackbody -I- powerlaw are (with 90% con- 
fidence error regions) nn = 0.22lg22 x 10 22 cm~ 2 , 
k,T = 0.19±g;g? keV. The implied 0.5-8 keV flux 
estimate is (2.8 ± 0.7) x 10~ 14 erg cm~ 2 s _1 with 
the error at the 68% confidence interval. This cor- 
responds to an unabsorbed flux of 4 x 10~ 14 erg 



of the source. The X-ray source properties are 
studied in detail in Camilo et al. (20091. The 



larger span of data included in this model results 
in significantly more counts in the light curve, con- 
firming the apparent 3-peak nature as proposed 
by Camilo et al. (20091, in contrast to the peak 



multiplicity of 2 assigned by Abdo et al. (2010d). 



3.9. PSR J1809 2332 

The timing model parameters for this pulsar 
are displayed in Table [12] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures 29] 30 and 31 respectively. 

This pulsar was discovered in the direction 
of the Galactic unidentified 7-ray source GeV 
J1809— 2327. Chandra observations revealed a 
probable pulsar with PWN that was proposed as 



the source of the 7-rays (|Braje et al. 2002|). The 
point source, 



CXOU J180950.2-233223, was as- 



sumed as the counterpart by Abdo et al. (2009a I 



As seen in Figure |29| the position error ellipse is 
very strongly elongated, again because of the very 
low ecliptic latitude of the source. Nevertheless, 
the X-ray point source is within the timing er- 
ror ellipse, strengthening the identification with 
this point source. If the TOAs are fitted with the 
position held fixed at the location of the Chan- 
dra point source, there are significant correlated 
residuals, which require a frequency third deriva- 
tive term in the model to give a reasonable x 2 f° r 
the fit. 



3.8. PSR J1741 2054 

The timing model parameters for this pulsar 
are displayed in Table [Tl] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [26] [27] and [28] respectively. 

The bright X-ray counterpart (Swift J174157.6— 
seen in Figure [26] was proposed as the likely coun- 
terpart to this pulsar in the discovery paper. Sub- 
sequently, a LAT timing position presented by 
Camilo et al. (2009), who also reported the dis- 



covery of radio pulsations from this pulsar, added 
confidence to this proposal, and the model we 
present here strengthens the case. The position 
error is still highly elongated in the declination 
direction because of the very low ecliptic latitude 



3.10. PSR J1813-1246 

The timing model parameters for this pulsar 
are displayed in Table [13] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures 32] 33 and 34 respectively. 
205411) This pulsar, which has the highest spindown 
luminosity of the first 16 blind search pulsars 
discovered with the LAT, exhibited a glitch on 
about 2009 September 20 with magnitude Au/u = 
1.17 x 10 -6 . The timing model for the glitch pre- 
sented here represents the data fairly well, but is 
not unique. Our model includes an instantaneous 
and permanent jump in the pulsar frequency and 
frequency derivative at the glitch. Other solutions 
are possible with slightly different glitch epochs 
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(within a day or so) and different parameters, or 
with decaying transient changes in the spin pa- 
rameters at the glitch. With 6 days between TOA 
measurements and significant timing noise seen in 
this pulsar, we are not able to distinguish between 
these possibilities. It is clear in Figure 33 that 



there are significant non-white residuals remain- 
ing in the data, indicating that the model is not 
fully accounting for the spindown behavior of the 
pulsar. With a longer span of post-glitch data a 
more definitive model may be possible. 

The bright X-ray source, Swift J181323.4-124600, 
was noted as the counterpart to this pulsar in the 
discovery paper and our timing confirms the asso- 
ciation. 

3.11. PSR J1826-1256 

The timing model parameters for this pulsar 



are displayed in Table 14 and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [35] [36] and [37] respectively. 

The fast spin and narrow pulse profile allow 
this pulsar to be localized to about 1". The tim- 
ing position is consistent with the X-ray point 
source AX J1826. 1-1257 (R.A. = 18:26:08.2, 
Decl. = —12:56:46), which was discovered in 
ASCA observations of the EGRET 7-ray source 
GeV J1825-1310 (Roberts et al.||2001b. An im- 



proved position of the X-ray point source (R.A. 
= 18:26:08.54, Decl. = -12:56:34.6) was derived 
from a Chandra image (M. Roberts, private com- 
munication), which is 1.6" from the timing posi- 
tion. The measured v for this pulsar is quite large, 
and the Monte Carlo error estimate yields 17", 
which is much larger than the offset seen between 
the timing position and the X-ray counterpart. 

3.12. PSR J1836+5925 

The timing model parameters for this pulsar 
are displayed in Table [15] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 



in Figures [38] [39] and [40] respectively. 

A detailed analysis and earlier timing model for 



this source were published by Abdo et al. (2010a I. 
Our model, which includes an additional 6 months 
of data is consistent with the earlier results with 
the addition of a weak detection of a frequency 



second derivative term. As seen in Figure [38] the 
timing position is fully consistent with the X-ray 
source RX J1836.2+5925, with the offset being 
only 0.2". 

3.13. PSR J1907+0602 

The timing model parameters for this pulsar 
are displayed in Table [16] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [44] [42] and [43] respectively. 

A detailed analysis of this pulsar, including an 
earlier timing solution and the discovery of radio 



pulsations, was presented by Abdo et al. (2010c I 



Our timing model is fully consistent within the 
errors to the one they presented, though we now 
have almost 5 months more data, which signif- 
icantly reduces the uncertainties in the parame- 
ters. The position reported for the Chandra point 
source is 2.3" from our best timing position, which 
is significantly larger than the 0.6" statistical er- 
ror in the timing position or the 0.6" error in the 
Chandra position. However, the offset is compara- 
ble to the expected systematic error from timing 
noise, based on the Monte Carlo simulations, so 
this is not strong evidence against the association. 

3.14. PSR J1958+2846 

The timing model parameters for this pulsar 
are displayed in Table [17] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [44] [45] and [46] respectively. 

In the discovery paper (Abdo et al. 2009a), the 
X-ray source Swift J195846.1+284602 was pro- 
posed as the likely counterpart and used for the 
name of the pulsar. As shown in Figure |44| the 
timing position no longer supports this identifica- 
tion, being offset from the X-ray source by 80". 
There is no significant X-ray counterpart detected 
at the timing position. There is also no indica- 
tion for strong timing noise in this pulsar, which 
might cause a large systematic error in the tim- 
ing position. Deeper X-ray observations will be 
required to detect the true X-ray counterpart for 
this source. 
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3.15. PSR J2021+4026 

The timing model parameters for this pulsar 
are displayed in Table [18] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures [47] [48] and [49] respectively. 

PSR J2021+4026 is the long-sought pulsar in 
the 7 Cygni supernova remnant (SNR G78. 2+2.1). 
In the discovery paper ( Abdo et al. [2009a ), it was 
pointed out that the X-ray source S21, as iden- 



tified earlier in Chandra observations (Weisskopf 
et al.|2006| |, was the most likely counterpart based 
on the initial pulsar timing. Our best fit position 
is 7.7" to the west of S21 (see Figure 471, an offset 
that is somewhat larger than the predicted sys- 
tematic error of 2.5" based on our Monte Carlo. 
However, when v is added to the model, the po- 
sition shifts by 4.6", so this is a lower bound on 
the systematic error from the timing noise. As 
S21 is still the closest X-ray source to the timing 
position, we conclude that it is indeed the likely 
counterpart. Longer term timing will improve our 
localization and reduce the systematic error con- 
tribution from timing noise. A similar conclusion 



was reached by Trepl et al. (2010). 



Both the timing position and the X-ray source 
S21 are well outside the 95% confidence localiza- 
tion ellipse of the LAT sources, as seen in Fig- 
ure [47] However, this region includes statistical 
errors only and this source is in the very compli- 
cated Cygnus region of the Galaxy. The local- 
ization of OFGL J2021. 5+4026 in the LAT Bright 



Source List (Abdo et al. 2009b) did include S21, 



but with the improved statistics using 18 months 
of data, systematic errors due to improperly mod- 
eled diffuse emission or unknown point sources can 
start to dominate the error budget. The associa- 
tion of the LAT source with the pulsar is beyond 
doubt because of the detection of pulsations. 

3.16. PSR J2032+4127 

The timing model parameters for this pulsar 
are displayed in Table [19] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 
in Figures 50] 51 and 52 respectively. 



This pulsar is studied in detail by |Camilo et al.| 
( 2009 ) , who reported the discovery of radio pulsa- 
tions from this source. The model presented here 



is consistent with theirs, and the positional asso- 
ciation with the Chandra point source MT91 213 



is confirmed (see Figure 50) 



3.17. PSR J2238+5904 

The timing model parameters for this pulsar 
are displayed in Table [20] and the timing position 
determination, post-fit residuals, 2-D phaseogram, 
and folded pulse profile for this pulsar are shown 



in Figures [53] [54] and [55] respectively. 

Even though this pulsar is quite faint, it can be 
timed with an RMS residual of 1 ms because of 
its very sharp pulse profile. Consequently, we now 
have a very precise timing position, but there is no 
significant X-ray counterpart detected in the Swift 



image (see Figure 53 1. The pulsar is located 0.6 



degree from the radio pulsar J2240+5832 detected 
recently in 7-rays (Theureau et al. 2010). The 



narrow pulse profile means that this pulsar can be 
blanked from the LAT data with a loss of only 
~ 20% of the exposure time. 

4. Radio Counterpart Searches 

All of these pulsars, except for PSR J1124-5916, 
were discovered in 7-ray searches and thus are 
7-ray-selected pulsars, but targeted radio obser- 
vations are required to determine if they are also 
radio quiet, or could have been discovered in radio 
surveys independently. The population statistics 
of radio-quiet vs. radio-loud 7-ray pulsars have 
important implications for 7-ray emission models 
(Gonthier et al. 2007). These observations are 



also important inputs into the population syn- 
thesis modeling of the full Galactic population 



of rotation-powered pulsars ( Faucher-Giguere & 
Kaspi 2006, for example). 



The precise positions derived from the LAT 
timing of these pulsars allowed us to perform deep 
follow up radio observations to search for pulsa- 
tions from each of the new pulsars. We used the 
NRAO 100-m Green Bank Telescope (GBT), the 
Arecibo 305-m radio telescope, and the Parkes 64- 
m radio telescope for these observations. The in- 
strument parameters used in the sensitivity calcu- 
lations are shown in Table[2] The log of obser- 
vations is shown in Table [3]and has columns for 
the target name, observation code (refer to Table 
[2]), observation date, observation duration (i bs), 
the R.A. and Decl. of the telescope pointing di- 
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rection, the offset from the true pulsar position, 
an estimate of the sky temperature in that direc- 
tion at the observing frequency, and our computed 
flux density limit (5 m i n ), as described below. The 
observations taken from the literature have S m i n 
recomputed in a consistent way as well as the orig- 
inally published flux limits in parentheses. In ad- 
dition, the fluxes of the detected pulsars are noted 
with "Det" in parentheses. It is notable that the 
flux of the detected pulsar J1907+0602 is below 
our nominal flux limit. This is caused primarily 
by the fact that the detected pulsar has a much 
smaller duty cycle than the 10% that we assume. 

All observations were taken in search mode 
(where all data are recorded without folding at a 
nominal pulse period) and the data were reduced 
using standard pulsar analysis software, such as 
PRESTO ( |Ransom et~aL||2002| >. In each case, we 
searched a range of dispersion measure (DM) trials 
out to a maximum of at least 2 times the maxi- 
mum DM value predicted by the NE2001 model 



(Cordes & Lazio 2002) for that direction. For 13 



of the 16 7-ray-selected pulsars, no radio pulsa- 
tions were detected and we report upper limits 
in Table [3] For three of the pulsars, pulsations 
were detected. Pulsations from PSRs J2032+4127 
and J1741— 2054 were reported by |Camilo et al 
(2009) and very faint radio pulsations from PSR 



J1907+06 were reported by Abdo et al. (2010c I 



Here, we compile the upper limits from the litera- 
ture as well as from our new observations. 

We calculate upper limits in a consistent man- 
ner for all of our observations as well as those from 
the literature. To compute the minimum pulsar 
flux that would have been detected in these obser- 
vations, we use the modified radiometer equation 
(e.g. Lorimer & Kramer 2005): 



5m in P 



w 



( S f N ) m i n T S y S 

GVVobxA/ V P-W 



(6) 



where j3 is the instrument-dependent factor due to 
digitization and other effects; (S/N) min = 5 is the 
threshold signal to noise for a pulsar to have been 
confidently detected; T sys = T rcc + T s ^ y , G is the 
telescope gain, n p is the number of polarizations 
used (2 in all cases); i bs is the integration time; 
A/ is the observation bandwidth; P is the pulsar 
period; W is the pulse width (for uniformity, we 
assume W = 0.1P). 

Because some of the observations were taken 



before the precise positions were known, some of 
the pointing directions are offset from the true di- 
rection to the pulsar. We use a simple approxi- 
mation of a telescope beam response to adjust the 
flux sensitivity in these cases . This factor is 



/ = cxp 



-(0/HWHM) 1 
L5 



(7) 



where is the offset from the beam center and 
HWHM is the beam half-width at half maximum. 
A computed flux limit of S at the beam center is 
thus corrected to S/f for a target offset from the 
pointing direction. The resultant flux limits are 
compiled in Table [3] 
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Table 2 

Definition of radio observing codes 



Obs Code 


Telescope 


Gain 


Freq 


A/ 


/3 a 


n p 


HWHM 


Tycc 






(K/Jy) 


(MHz) 


(MHz) 






(arcmin) 


(K) 


GBT-350 


GBT 


2.0 


350 


100 


1.05 


2 


18.5 


46 


GBT-820 


GBT 


2.0 


820 


200 


1.05 


2 


7.9 


29 


GBT-820BCPM 


GBT 


2.0 


820 


48 


1.05 


2 


7.9 


29 


GBT-S 


GBT 


1.9 


2000 


700 b 


1.05 


2 


3.1 


22 


AO-327 


Arecibo 


11 


327 


50 


1.12 


2 


6.3 


116 


AO-430 


Arecibo 


11 


430 


40 


1.12 


2 


4.8 


84 


AO-Lwide 


Arecibo 


10 


1510 


300 


1.12 


2 


1.5 


27 


AO- ALFA 


Arecibo 


10 


1400 


100 


1.12 


2 


1.5 


30 


Parkes-MB256 


Parkes 


0.735 


1390 


256 


1.25 


2 


7.0 


25 


Parkes-AFB 


Parkes 


0.735 


1374 


288 


1.25 


2 


7.0 


25 


Parkes-BPSR 


Parkes 


0.735 


1352 


340 


1.05 


2 


7.0 


25 



a Instrument-dependent sensitivity degradation factor, see equation [6] 

b The instrument records 800 MHz of bandwidth, but to account for a notch filter for RFI 
and the lower sensitivity near the band edges, we use an effective bandwidth of 700 MHz for 
the sensitivity calculations. 
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Fig. 5.— Timing position for PSR J0007+7303 
(blue ellipse). The large green ellipse is 
the LAT 95% confidence localization of 1FGL 
J0007.0+7303, based on 18 months of data. 
The red X is the Chandra position of RX 
J0007. 0+ 7303, which is 4.4" f rom the best timing 
position ( Halpern et al.|[2004 ). 
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Fig. 7. — 2-D phaseogram and pulse profile of PSR 
J0007+7303. Two rotations are shown on the X- 
axis. The photons were selected according to the 
ROI and E m i n in Table |4j In this and the follow- 
ing phaseogram plots, the grey scale is the num- 
ber of photons in each phase/time bin, without 
any correction for exposure, so apparent bright- 
ness changes are caused by the precession period 
of the Fermi satellite, interruptions in science op- 
erations, or from operational changes in the rock- 
ing pattern in sky survey mode. The fiducial point 
corresponding to TZRMJD is phase 0.0. This and 
all following pulse profiles arc constructed with 32 
bins across the pulse period. 
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Fig. 6. — Post-fit timing residuals for PSR 
J0007+7303. 
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Table 4 
PSR J0007+7303 



Parameter 




Value 


Right ascension, a (J2000.0) . . 




00:07:00.6 ±0.2 S 


Declination, S (J2000.0) 




+73:03:07.0 ±0.6" 


Monte Carlo position uncertainty 


2" 


Pulse frequency, v (s _1 ) 




3.165827380(3) 


Frequency first derivative, v (s 


- 2 )-. 


-3.6136(2)xl0- 12 


Frequency second derivative, v 


(s- 3 ) 


-7(1) xlO- 23 


Epoch of Frequency (MJD) 




54952 


Glitch Epoch 




54952.652 


Glitch Av (s' 1 ) 




1.759(3) xlO" 6 


Glitch Av (s- 2 ) 







TZRMJD 




54952.334185720257651 


Number of photons (n 7 ) 




12790 


Number of TOAs 




55 


RMS timing residual (ms) 




2.2 


Template Profile 




KDE 


K ■ 




150 McV 


ROI 




1.5° 


Valid range (MJD) 




54682 - 55222 
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Fig. 8.— Timing position for PSR J0357+3205 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J0357.8+3205, based on 18 
months of data. The background 0.2-10 keV X- 
ray image is from a 2.6 ks Swift observation (Ob- 
sID 00031299001), smoothed with a gaussian with 
o = 7". 
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Fig. 10. — 2-D phaseogram and pulse profile of 
PSR J0357+3205. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E m i n in Table[5j The fiducial point 
corresponding to TZRMJD is phase 0.0. 
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Fig. 9. — Post-fit timing residuals for PSR 
J0357+3205. 
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Table 5 
PSR J0357+3205 



Parameter 


Value 


Right ascension, a (J2000.0) 


03:57:52.5 ±0.2 S 


Declination, 5 (J2000.0) 


+32:05:25 ±6" 


Monte Carlo position uncertainty 


18" 


Pulse frequency, v (s~ x ) 


2.251722292(3) 


Frequency first derivative, v (s~ 2 ) . . 


-6.61(l)xl0~ 14 


Frequency second derivative, v (s~ 3 ) 


\i>\ < 6 x 10~ 23 


Epoch of Frequency (MJD) 


54946 


TZRMJD 


54946.341346723796502 


Number of photons (n 7 ) 


1335 


Number of TOAs 


25 


RMS timing residual (ms) 


5.3 


Template Profile 


1 Gaussian 


E ■ 


250 McV 


ROI 


0.8° 


Valid range (MJD) 


54682 - 55210 
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Fig. 11— Timing position for PSR J0633+0632 
(blue ellipse). The large ellipse is the LAT posi- 
tion of 1FGL J0633.7+0632, based on 18 months 
of data. The background 0.5-8 keV X-ray image 
is a 20 ks Chandra ACIS-S image (ObsID 11123), 
smoothed with a gaussian with a = 0.5". The in- 
set shows a 10" region around the timing location. 
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Pulse Phase 

Fig. 13. — 2-D phaseogram and pulse profile of 
PSR J0633+0632. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E m i n in Table[6j The fiducial point 
corresponding to TZRMJD is phase 0.0. 



Fig. 12. — Post-fit timing residuals for PSR 
J0633+0632. 



24 



Table 6 
PSR J0633+0632 



Parameter 


Value 


Right ascension, a (J2000.0) 


06:33:44.21 ±0.02 s 


Declination, 5 (J2000.0) 


+06:32:34.9 ±1.6" 


Monte Carlo position uncertainty 


3.5" 


Pulse frequency, v (s~ x ) 


3.3625291588(7) 


Frequency first derivative, v (s~ 2 ) . . 


-8.9991(3)xl0- 13 


Frequency second derivative, v (s~ 3 ) 


-2(l)xl0- 23 


Epoch of Frequency (MJD) 


54945 


TZRMJD 


54945.385967311181439 


Number of photons (n 7 ) 


1174 


Number of TOAs 


23 


RMS timing residual (ms) 


1.4 


Template Profile 


2 Gaussian 


E ■ 


550 McV 


ROI 


0.6° 


Valid range (MJD) 


54682 - 55208 
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Fig. 14.— Timing position for PSR J1124-5916 
(small blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1124.6-5916, based on 18 
months of data. The red cross marks the position 



of the Chandra point source ( Camilo et al. 2002 ) 



associated with the pulsar (see [ 3.4 1 
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Fig. 15. — Post-fit timing residuals for PSR J1124- 
5916. 



Fig. 16. — 2-D phaseogram and pulse profile of 
PSR J1124— 5916. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E m i n in Table[7J The fiducial point 
corresponding to TZRMJD is phase 0.0. The red 
line is a 1.4 GHz radio profile from the Parkes 
radio telescope, with the correct absolute phase 
alignment. 
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Table 7 
PSR J1124-5916 



Parameter Value 



Right ascension, a (J2000.0) 11:24:39.0(1) 

Declination, S (J2000.0) -59:16:19(1) 

Pulse frequency, v (s^ 1 ) 7.381334652(9) 

Frequency first derivative, i> (s~ 2 ) . . -4.10029(9) x 10 -11 

Frequency second derivative, v (s~ 3 ) — 8.6(4) x 10~ 22 

Epoch of Frequency (MJD) 54683.281414 

Glitch Epoch (MJD) 55191 

Glitch Av (s- 1 ) 1.18(9) xl0~ 7 

Glitch Ai> (s~ 2 ) 1.94(2) xl0~ 13 

TZRMJD 55053.0521054597626 

TZRFREQ (MHz) 1371.067 

TZRSITE 7 (Parkes) 

Number of photons (n 7 ) 5030 

Number of TOAs 40 

RMS timing residual (ms) 2.8 

Template Profile 2 Gaussian 

E min 200 MeV 

ROI 0.9° 

Valid range 54682-55415 
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Fig. 17.— Timing position for PSR J1418-6058 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1418.7-6057, based on 
18 months of data. Red crosses mark the posi- 
tions of the sources Rl and R2 (see ^3.5). The 
X-ray image is a 70 ks Chandra ACIS observation 
fro m 2007 July 14 (ObsID 7640), first published 
by Roberts ( 2008[ ) before the pulsar itself was de- 
tected. 
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Fig. 19. — 2-D phaseogram and pulse profile of 
PSR J1418-6058. Two rotations are shown on the 
X-axis. The photons were selected according to 
the ROI and E m i n in Table [§} The fiducial point 
corresponding to TZRMJD is phase 0.0. 



Fig. 18.— Post-fit timing residuals for PSR J1418- 
6058. 
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Table 8 
PSR J1418-6058 



Parameter 


Value 


Right ascension, a (J2000.0) 


f4:f8:42.7 ±0.1" 


Declination, <5 (J2000.0) 


-60:57:49 ±2" 


Monte Carlo position uncertainty 


7" 


Pulse frequency, v (s _1 ) 


9.043798163(1) 


Frequency first derivative, i> (s~ 2 ) . . 


-1.38548(8) xlO- 11 


Frequency second derivative, v (s~ 3 ) 


6.4(3) xl0~ 22 


Frequency third derivative, v (s~ 4 ) . 


-8(2) xlO- 29 


Epoch of Frequency (MJD) 


54944 


TZRMJD 


54944.2886329214 


Number of photons (n 7 ) 


7283 


Number of TOAs 


33 


RMS timing residual (ms) 


1.9 


Template Profile 


KDE 


E ■ 


250 MeV 


ROl 


0.5° 


Valid range (MJD) 


54682 - 55205 
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Fig. 20.— Timing position for PSR J1459-6053 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1459.4-6053, based on 
18 months of data. The background 0.2-10 keV 
X-ray image is a 6.8 ks Swift image (ObsID 
00031359002), smoothed with a gaussian with a = 
7" . A 10 x 10" region around the timing position 
is indicated with the black square and shown in 
the inset at the lower left. 
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Fig. 22. — 2-D phaseogram and pulse profile of 
PSR J1459— 6053. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and i^mm in Table[9j The hducial point 
corresponding to TZRMJD is phase 0.0. 



Fig. 21.— Post-fit timing residuals for PSR J1459- 
6053. 
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Table 9 
PSR J1459-6053 



Parameter 


Value 


Right ascension, a (J2000.0) 


14:59:29.99 ±0.06 s 


Declination, <5 (J2000.0) 


-60:53:20.7 ±0.4" 


Monte Carlo position uncertainty 


1.3" 


Pulse frequency, v (s _1 ) 


9.694559498(1) 


Frequency first derivative, v (s~ 2 ) . . 


-2.37503(5)xl0~ 12 


Frequency second derivative, v (s~ 3 ) 


-4(2) xl0~ 23 


Epoch of Frequency (MJD) 


54935 


TZRMJD 


54936.19962194 


Number of photons (n 7 ) 


3305 


Number of TOAs 


26 


RMS timing residual (ms) 


1.1 


Template Profile 


KDE 


R ■ 


350 McV 


ROI 


0.7° 


Valid range (MJD) 


54682-55210 
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Fig. 23.— Timing position for PSR J1732-3131 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1732. 5-3131, based on 18 
months of data. The background 0.5-8 keV X-ray 
image is a 20 ks Chandra ACIS-S image (ObsID 
11125), smoothed with a gaussian with a = 0.5". 




1.0 

Pulse Phase 



Fig. 25. — 2-D phaseogram and pulse profile of 
PSR J1732— 3131. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [lOj The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 24. — Post-fit timing residuals for PSR 
J1732-3131. 
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Table 10 
PSR J1732-3131 



Parameter Value 



Right ascension, a (J2000.0) 17:32:33.54 ±0.03' 

Declination, <5 (J2000.0) -31:31:23 ±2" 

Monte Carlo position uncertainty 3" 

Pulse frequency, v (s^ 1 ) 5.0879411200(5) 

Frequency first derivative, v (s~ 2 ) . . — 7.2609(3)x \{)~ Vi 

Frequency second derivative, v (s -3 ) \v\ < 2 x 10~ 23 

Epoch of Frequency (MJD) 54933 

TZRMJD 54957.3282196892 

Number of photons (n 7 ) 4236 

Number of TOAs 22 

RMS timing residual (ms) 1.0 

Template Profile KDE 

E min 400 McV 

ROI 0.5° 

Valid range (MJD) 54682-55207 
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Fig. 26.— Timing position for PSR J1741-2054 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1741.8-2101, based on 
18 months of data. The background 0.2-10 keV 
X-ray image is a 4.3 ks Swift XRT image (Ob- 
sID 00031277001), smoothed with a gaussian with 
o = 7". 
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Fig. 28. — 2-D phaseogram and pulse profile of 
PSR J1741-2054. Two rotations are shown on the 
X-axis. The photons were selected according to 
the ROI and Em\n 



in Table 11 The fiducial point 
corresponding to TZRMJD is phase 0.0. 
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Fig. 27. — Post-fit timing residuals for PSR 
J1741-2054. 
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Table 11 
PSR J1741-2054 



Parameter Value 



Right ascension, a (J2000.0) 17:41:57.23 ±0.05' 

Declination, 5 (J2000.0) -20:53:57 ±19" 

Monte Carlo position uncertainty 20" 

Pulse frequency, v (s^ 1 ) 2.417209833(1) 

Frequency first derivative, v (s~ 2 ) . . —9.923(3) x 10 -14 

Frequency second derivative, v (s~ 3 ) \v\ < 2 x 10 -23 

Epoch of Frequency (MJD) 54933 

TZRMJD 54945.3859666189 

Number of photons (n 7 ) 3135 

Number of TOAs 23 

RMS timing residual (ms) 2.6 

Template Profile KDE 

E min 300 McV 

ROI 0.8° 

Valid range (MJD) 54682-55208 
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Fig. 29.— Timing position for PSR J1809-2332 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1809.8-2332, based on 18 
months of data. The background image is a 9.8 ks 
Chandra ACIS-I image (Obsld 739) showing the 
bright point source CXOU J180950. 2-233223. 
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Fig. 31. — 2-D phaseogram and pulse profile of 
PSR J1809— 2332. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [12] The hducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 30. — Post-fit timing residuals for PSR 
J1809-2332. 
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Table 12 
PSR J1809-2332 



Parameter 


Value 


Right ascension, a (J2000.0) 


18:09:50.31 ±0.06 s 


Declination, <5 (J2000.0) 


-23:33:35 ±51" 


Monte Carlo position uncertainty 


28" 


Pulse frequency, v (s _1 ) 


6.8125205463(3) 


Frequency first derivative, v (s~ 2 ) . . 


-1.59748(l)xl0~ 12 


Frequency second derivative, v (s~ 3 ) 


\i>\ < 1 x 10~ 23 


Epoch of Frequency (MJD) 


54935 


TZRMJD 


54947.1551911038 


Number of photons (n 7 ) 


10422 


Number of TOAs 


27 


RMS timing residual (ms) 


0.4 


Template Profile 


KDE 


E ■ 


250 McV 


ROI 


0.8° 


Valid range (MJD) 


54682-55211 
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Fig. 32.— Timing position for PSR J1813-1246 
(yellow ellipse). The large green ellipse is the 
LAT position of 1FGL J1813.3-1246, based on 
18 months of data. The background X-ray im- 
age is a 3.2 ks Swift XRT observation (Ob- 
sID 00031381001), smoothed with a gaussian of 
width 7", showing the bright point source Swift 
J181323.4-124600. 
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Fig. 34. — 2-D phaseogram and pulse profile of 
PSR J1813— 1246. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [13] The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 33. — Post-fit timing residuals for PSR 
J1813-1246. 
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Table 13 
PSR J1813-1246 



Parameter Value 



Right ascension, a (J2000.0) 18:13:23.77 ±0.01 

Declination, <5 (J2000.0) -12:45:59.2 ±1.5' 

Monte Carlo position uncertainty 1.5" 

Pulse frequency, v (s^ 1 ) 20.802023359(3) 

Frequency first derivative, v (s~ 2 ) . . — 7.60023(9)x 10~ 

Frequency second derivative, v (s~ 3 ) \v\ < 6 x 10~ 23 

Glitch Epoch 55094.1227 

Glitch Av (s- 1 ) 2.4256(9) xlO" 5 

Glitch Az> (s- 2 ) -4.9(2) x 10" 14 

TZRMJD 54954.309848738 

Number of photons (n 7 ) 11611 

Number of TOAs 91 

RMS timing residual (ms) 0.7 

Template Profile KDE 

E min 200 McV 

ROI 1.0° 

Valid range 54682-55226 
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Fig. 35.— Timing position for PSR J1826-1256 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1826.1-1256, based on 
18 months of data. The background 0.2-10 keV 
X-ray image is a 4.3 ks Swift XRT image (Ob- 
sID 00035179002), smoothed with a gaussian with 
a = 7". The red cross marks the Chandra position 
of X-ray point source AX J1826.1-1257. 
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Fig. 37. — 2-D phaseogram and pulse profile of 
PSR J1826— 1256. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [l4j The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 36. — Post-fit timing residuals for PSR 
J1826-1256 
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Table 14 
PSR J1826-1256 



Parameter 


Value 


Right ascension, a (J2000.0) 


18:26:08.53 ±0.01 s 


Declination, <5 (J2000.0) 


-12:56:33.0 ±0.5" 


Monte Carlo position uncertainty 


17" 


Pulse frequency, v (s _1 ) 


9.0724588059(3) 


Frequency first derivative, v (s~ 2 ) . . 


-9.99654(l)xl0~ 12 


Frequency second derivative, v (s~ 3 ) 


1.85(5) xl0~ 22 


Epoch of Frequency (MJD) 


54934 


TZRMJD 


54946.3413482956 


Number of photons (n 7 ) 


10860 


Number of TOAs 


25 


RMS timing residual (ms) 


0.28 


Template Profile 


KDE 


R ■ 


200 McV 


ROI 


0.6° 


Valid range 


54682-55210 
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Fig. 38.— Timing position for PSR J1836+5925 
(yellow ellipse). The large green ellipse is the 
LAT position of 1FGL J1836.2+5925, based on 
18 months of data. The X-ray image is a 46 ks 
Chandra HRC image (Obsld 6182) and the point 
source at the timing position is RX J1836. 2+5925. 
The inset (3.0" in width) shows the region of the 
source in more detail. 
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Fig. 40. — 2-D phaseogram and pulse profile of 
PSR J1836+5925. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [15] The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 39. — Post-fit timing residuals for PSR 
J1836+5925. 
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Table 15 
PSR J1836+5925 



Parameter Value 



Right ascension, a (J2000.0) 18:36:13.69 ±0.02 s 

Declination, 5 (J2000.0) +59:25:30.0 ±0.3" 

Monte Carlo position uncertainty < 3" 

Pulse frequency, v (s" 1 ) 5.7715509149(4) 

Frequency first derivative, v (s~ 2 ) . . —5.007(2) x 10 -14 

Frequency second derivative, v (s~ 3 ) — 1.3(8) xl0~ 23 

Epoch of Frequency (MJD) 54935 

TZRMJD 54936.1996219705 

Number of photons (n 7 ) 14875 

Number of TOAs 26 

RMS timing residual (ms) 0.85 

Template Profile KDE 

E min 200 McV 

ROI 1.6° 

Valid range 54682-55210 
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Fig. 41.— Timing position for PSR J1907+0602 
(yellow ellipse). The large green ellipse is the 
LAT position of 1FGL J1907.9+0602, based on 
18 months of data. The X-ray image is a 0.75- 
2.0 kcV Chandra image (ObsID 11124). The inset 
(7.2" in width) shows a detail of the region around 
the pulsar. 
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Fig. 43. — 2-D phaseogram and pulse profile of 
PSR J1907+0602. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [16] The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 42. — Post-fit timing residuals for PSR 
J1907+0602. 
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Table 16 
PSR J1907+0602 



Parameter Value 



Right ascension, a (J2000.0) 19:07:54.74 ±0.01 s 

Declination, 5 (J2000.0) +06:02:16.9 ±0.3" 

Monte Carlo position uncertainty 2.5" 

Pulse frequency, v (s^ 1 ) 9.3779822336(4) 

Frequency first derivative, v (s~ 2 ) . . -7.63559(2) x 10 -12 

Frequency second derivative, v (s~ 3 ) 1.88(7) x 10 -22 

Epoch of Frequency (MJD) 54935 

TZRMJD 54947.1551911789 

Number of photons (n 7 ) 10629 

Number of TOAs 27 

RMS timing residual (ms) 0.47 

Template Profile KDE 

E min 50 MeV 

ROI 0.7° 

Valid range 54682-55211 
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Fig. 44.— Timing position for PSR J1958+2846 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J1958.6+2845, based on 
18 months of data. The background 0.2-10 keV 
X-ray image is a 5.5 ks Swift XRT image (Ob- 
sID 00031374001), smoothed with a gaussian with 
a = 7". The red 'x' marks the position of Swift 
J195846.1+284602. 
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Fig. 46. — 2-D phaseogram and pulse profile of 
PSR J1958+2846. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table (T7] The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 45. — Post-fit timing residuals for PSR 
J1958+2846. 
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Table 17 
PSR J1958+2846 



Parameter Value 



Right ascension, a (J2000.0) 19:58:40.07 ±0.03' 

Declination, 5 (J2000.0) +28:45:54 ±1" 

Monte Carlo position uncertainty 3.5" 

Pulse frequency, v (s" 1 ) 3.4436537099(5) 

Frequency first derivative, v (s~ 2 ) . . -2.5145(2) x 10 -12 

Frequency second derivative, v (s~ 3 ) 3(2) xlO -23 

Epoch of Frequency (MJD) 54800 

TZRMJD 54957.3282188686 

Number of photons (n 7 ) 1910 

Number of TOAs 26 

RMS timing residual (ms) 2.1 

Template Profile 2 Gaussian 

E min 550 McV 

ROI 0.6° 

Valid range 54682-55210 
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Fig. 47.— Timing position for PSR J2021+4026 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J2021.5+4026, based on 
18 months of data. The red 'x' marks the position 
of the source S21 (see text). The background X- 
ray image is a portion of a Chandra ACIS-I image 
(ObsID 5533), with 3-pixel gaussian smoothing. 
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Fig. 49. — 2-D phaseogram and pulse profile of 
PSR J202 1+4026. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [18] The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 48. — Post-fit timing residuals for PSR 
J2021+4026. 
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Table 18 
PSR J2021+4026 



Parameter Value 



Right ascension, a (J2000.0) 20:21:29.99 ±0.03 s 

Declination, 5 (J2000.0) +40:26:45.1 ±0.7" 

Monte Carlo position uncertainty 2.5" 

Pulse frequency, v (s" 1 ) 3.7690668480(6) 

Frequency first derivative, v (s~ 2 ) . . -7.7681(3) x 10 -13 

Frequency second derivative, v (s~ 3 ) 3.9(2) xl0~ 22 

Epoch of Frequency (MJD) 54936 

TZRMJD 54957.3282196715 

Number of photons (n 7 ) 11853 

Number of TOAs 30 

RMS timing residual (ms) 2.0 

Template Profile KDE 

E min 400 McV 

ROI 0.7° 

Valid range 54682-55213 
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Fig. 50.— Timing position for PSR J2032+4127 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J2032.2+4127, based on 
18 months of data. The background image is from 
a 49 ks Chandra ACIS observation (ObsID 4501). 
The inset shows a 3.6" region around the timing 
position in more detail. 
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Fig. 52. — 2-D phaseogram and pulse profile of 
PSR J2032+4127. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [19] The fiducial 
point corresponding to TZRMJD is phase 0.0. 



Fig. 51. — Post-fit timing residuals for PSR 
J2032+4127. 
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Table 19 
PSR J2032+4127 



Parameter Value 



Right ascension, a (J2000.0) 20:32:13.25 ±0.01 s 

Declination, 5 (J2000.0) +41:27:24.8 ±0.3" 

Monte Carlo position uncertainty 3" 

Pulse frequency, v (s^ 1 ) 6.9809196293(4) 

Frequency first derivative, v (s~ 2 ) . . -9.9293(2)x 10~ 13 

Frequency second derivative, v (s~ 3 ) —1.88(1) x 10 -21 

Epoch of Frequency (MJD) 54938 

TZRMJD 54951.224402859 

Number of photons (n 7 ) 1633 

Number of TOAs 45 

RMS timing residual (ms) 0.9 

Template Profile 2 Gaussian 

E min 900 McV 

ROI 0.5° 

Valid range 54682-55220 
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Fig. 53.— Timing position for PSR J2238+5903 
(blue ellipse). The large green ellipse is the 
LAT position of 1FGL J2238.4+5903, based on 
18 months of data. The background 0.2-10 keV 
X-ray image is a 4.9 ks Swift image (ObsID 
00031398001), smoothed with a gaussian with a = 
7". The inset shows a 7.2" region around the pul- 
sar in more detail. 
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Fig. 55. — 2-D phaseogram and pulse profile of 
PSR J2238+5903. Two rotations are shown on 
the X-axis. The photons were selected according 
to the ROI and E min in Table [20] The fiducial 
point corresponding to TZRMJD is phase 0.0. 
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Fig. 54. — Post-fit timing residuals for PSR 
J2238+5903. 
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Table 20 
PSR J2238+5903 



Parameter Value 



Right ascension, a (J2000.0) 22:38:28.27 ±0.04 s 

Declination, 5 (J2000.0) +59:03:40.8 ±0.4" 

Monte Carlo position uncertainty 3" 

Pulse frequency, v (s^ 1 ) 6.1450029089(4) 

Frequency first derivative, v (s~ 2 ) . . — 3.6641(2)x 10~ 12 

Frequency second derivative, v (s~ 3 ) 1.1(2) xlO -22 

Epoch of Frequency (MJD) 54800 

TZRMJD 54947.1551907197 

Number of photons (n 7 ) 1697 

Number of TOAs 27 

RMS timing residual (/is) 1171 

Template Profile 2 Gaussian 

E min 250 McV 

ROI 0.5° 

Valid range 54682-55211 
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5. Discussion 

We have presented precise phase-coherent tim- 
ing models using 18 months of Fermi LAT data for 
17 radio-quiet or radio- faint 7-ray pulsars. This 
provides precise spin parameters for all of the pul- 
sars and position determinations of order arcsec- 
ond accuracy. 

In most cases the new position determinations 
served to confirm previously proposed X-ray coun- 
terparts. However, in one case (PSR J1958+2846) 
the previously proposed counterpart is strongly 
discrepant with the new position determination. 
In one other case (PSR J1459— 6053), an X-ray 
source is apparent in a Swift image at the pulsar 
position. Lastly in 3 cases (PSRs J0633+0632, 
J1418-6058, and J2021+4026), the situation is a 
bit more complicated, because the observed off- 
sets between the timing position and the X-ray 
counterpart position may be accounted for by the 
effects of timing noise on the model fits. These 
were covered on a case-by-case basis. 

In three of the 17 pulsars (PSRs J0007+7303, 
J1124-5916, and J1813-1246), we have detected 
a glitch. This is not unexpected for a population 
of mostly young pulsars with characteristic ages 
of 10 4 - 10 5 years. These three glitches observed 
in the seventeen radio quiet and radio faint pul- 
sars are typical of the eight glitches observed in 
eighteen months of the Fermi 7-ray pulsars. All 
pulsars observed to glitch with the LAT pulsars 
have spin down energies E > 4.5 x 10 35 erg s _1 . In 
fact, most (6/8) of the glitching pulsars are above 
E>lx 10 36 erg s" 1 . All of the LAT glitching pul- 
sars have characteristic ages between 1-100 kyr. A 
more detailed analysis of timing across glitches in 
7-ray pulsars is in preparation (Dormody et al. 
2010, in preparation). 

With 18 months of timing data, we also have 
measurements of v for most of the pulsars. The 
measured vs are dominated by timing noise rather 
than the secular spin down behavior of the pulsars. 



Previously, Arzoumanian et al. (1994b) have de- 



fined a pulsar stability parameter 

1 

6~^' 



A W= lo giohd^ 



(8) 



where t is the observation duration and they define 
A 8 = A(10 8 s). They find a correlation of this 
stability parameter with pulsar period derivative 



(P), with the form 

A 8 =6.6 + 0.61og 10 P. (9) 
This relationship has been re-fit using a larger 



sample of pulsars by Hobbs et al. (2010), who ob- 



tain the following parameters: 

A 8 = 5.1 + 0.51og 10 P. (10) 
s of data, so we compute our 



10 s. As seen in Figure 56 



We do not have 10 
A parameter at t = 
we see a similar correlation with period derivative, 
albeit with a large amount of scatter. The two 
pulsars that stand out farthest from the relation 
as having very large v for their period derivatives 
are PSRs J2021+4026 and J2032+4127. 
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Fig. 56. — A parameter characterizing timing 
noise vs. v for these pulsars. The red arrows rep- 
resent 2-tr upper limits. The solid green line is 
the relation for A 8 found by | Arzoumanian et al.| 



( 1994b ) , while the dashed line is the relation found 



by Hobbs et al. (2010). Note that our data span 
is about half of the 10 8 s used in the definition of 
A 8 . 

These pulsars will continue to be timed regu- 
larly throughout the LAT mission. Using the fake 
TOA simulation capability of Tempo2, we have 
evaluated the possibility of measuring further as- 
trometric parameters for these pulsars. We find 
that in a 10 year mission, we are unlikely to be 
able to detect parallax or proper motion for any 
of these sources. The most nearby pulsar is ~ 400 
pc distant, and a parallax signal at that distance 
is 3 /zs, for a pulsar at an ecliptic latitude of 0. Un- 
fortunately, the nearby pulsars are also those with 
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the longest periods and the largest RMS timing 
residuals. We evaluated the possibility of detect- 
ing proper motion for a large transverse velocity 
of 1000 km s -1 , and again found that none of the 
pulsars look like promising candidates for proper 
motion measurements within the Fermi mission. 

We also made deep searches for radio pulsa- 
tions from the 7-ray selected pulsars. We compare 
these flux limits with the measured fluxes of the 
population of pulsars in the ATNF pulsar catalog 
(Manchester et al. 2005) in Figure 57 To make 



the fluxes comparable, we have scaled them all to 
the equivalent 1400 MHz flux density using a typ- 
ical pulsar spectral index of 1.6. The upper limits 
we have obtained are comparable to some of the 
faintest known radio pulsars, but the discovery of 
3.5 /xJy puls ations from PSR J1907+0602 ( [ArJdo 
et al. 2010c) raises the possibility that some of 



these could yet be detected in even deeper radio 
searches. 

The radio upper limits for 8 new 7-ray selected 



pulsar discovered with Fermi are presented in |Saz 
Parkinson et al. (2010). When combined with 
the results presented here, we now have deep up- 
per limits on all known 7-ray selected pulsars. 
A discussion of the radio upper limits on PSR 



J1836+5925 was also presented by |Abdo et al. 
pHOa 
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Fig. 57. — Summary of radio pulsation searches of 
7-ray selected pulsars. The red arrows denote up- 
per limits for the non-detections reported in Ta- 
ble [3] The large green dots are the radio pul- 
sation detections that have been previously re- 
ported ( jCamilo et al 



2009 Abdo et al. 2010c) 



The black dots are the 1400 MHz flux densities of 
the non-millisecond pulsars in the ATNF catalog, 
for comparison. The detections and upper limits 
that were made at different observing frequencies 
were scaled to 1400 MHz using a typical pulsar 
spectral index of 1.6. 



55 



(JAXA) in Japan, and the K. A. Wallenberg Foun- 
dation, the Swedish Research Council and the 
Swedish National Space Board in Sweden. 

Additional support for science analysis during 
the operations phase is gratefully acknowledged 
from the Istituto Nazionale di Astrofisica in Italy 
and the Centre National d'Etudes Spatiales in 
France. 

The Arecibo Observatory is part of the National 
Astronomy and Ionosphere Center, which is op- 
erated by Cornell University under a cooperative 
agreement with the National Science Foundation. 
The National Radio Astronomy Observatory is a 
facility of the National Science Foundation oper- 
ated under cooperative agreement by Associated 
Universities, Inc. The Parkes Observatory is part 
of the Australia Telescope which is funded by the 
Commonwealth of Australia for operation as a Na- 
tional Facility managed by CSIRO. 

Facilities: Fermi LAT, GBT, Arecibo, Parkes 



•50 



REFERENCES 

Abdo, A. A., et al. 2008, Science, 322, 1218, 
(CTA1) 

— . 2009a, Science, 325, 840, (16 Blind Search Pul- 
sars) 

— . 2009b, ApJS, 183, 46, (LAT Bright Source 
List) 

— . 2009c, Science, 326, 1512, (Cyg X-3) 

— . 2010a, ApJ, 712, 1209, (PSR J1836+5925) 

— . 2010b, ApJ, 720, 272, (Geminga) 

— . 2010c, ApJ, 711, 64, (PSR J1907+0602) 

— . 2010d, ApJS, 187, 460, (LAT Pulsar Catalog) 

— . 2010e, ApJ, 713, 154, (Vela2) 

— . 2011, ApJS, in press (Erratum to ApJS, 187, 
460) 

Andersson, N., Comer, G. L., & Prix, R. 2003, 
Physical Review Letters, 90, 091101 

Arzoumanian, Z., Nice, D. J., Taylor, J. FL, & 
Thorsett, S. E. 1994a, ApJ, 422, 671 

— . 1994b, ApJ, 422, 671 

Atwood, W. B., et al. 2009, ApJ, 697, 1071 

Becker, W., et al. 2004, ApJ, 615, 897 

Braje, T. M., Romani, R. W., Roberts, M. S. E., 
& Kawai, N. 2002, ApJ, 565, L91 

Camilo, F., Manchester, R. N., Gaensler, B. M., 
Lorimer, D. R., & Sarkissian, J. 2002, ApJ, 567, 
L71 

Camilo, F., et al. 2009, ApJ, 705, 1 

Cordes, J. M., & Downs, G. S. 1985, ApJS, 59, 
343 

Cordes, J. M., & Lazio, T. J. W. 2002, arXiv:astro- 
ph/0207156 

de Jager, O. C, Raubenheimer, B. C, & 
Swanepoel, J. W. H. 1986, A&A, 170, 187 

— . 1989, A&A, 221, 180 



Edwards, R. T., Hobbs, G. B., & Manchester, 
R. N. 2006, MNRAS, 372, 1549 

Faucher-Giguere, C, & Kaspi, V. M. 2006, ApJ, 
643, 332 

Gonthier, P. L., Story, S. A., Clow, B. D., & Hard- 
ing, A. K. 2007, Ap&SS, 309, 245 

Halpern, J. P., Camilo, F., & Gotthelf, E. V. 2007, 
ApJ, 668, 1154 

Halpern, J. P., Gotthelf, E. V., Camilo, F., 
Helfand, D. J., & Ransom, S. M. 2004, ApJ, 
612, 398 

Hobbs, G., Lyne, A. G., & Kramer, M. 2010, MN- 
RAS, 402, 1027 

Hobbs, G., Lyne, A. G., Kramer, M., Martin, 
C. E., & Jordan, C. 2004, MNRAS, 353, 1311 

Hobbs, G. B., Edwards, R. T., & Manchester, 
R. N. 2006, MNRAS, 369, 655 

Jackson, M. S., & Halpern, J. P. 2005, ApJ, 633, 
1114 

Kramer, M., & Wex, N. 2009, Classical and Quan- 
tum Gravity, 26, 073001 

Link, B., Epstein, R. I., & Lattimer, J. M. 1999, 
Physical Review Letters, 83, 3362 

Livingstone, M. A., Ransom, S., Camilo, F., 
Kaspi, V. M., Lyne, A., Kramer, M., & 
Stairs, I. H. 2009, ApJ, 706, 1163, submitted 
(arXiv:0901.2119) 

Lorimer, D., & Kramer, M. 2005, Handbook 
of Pulsar Astronomy (Cambridge University 
Press) 

Manchester, R. N., Hobbs, G. B., Teoh, A., & 
Hobbs, M. 2005, AJ, 129, 1993 

Ng, C, Roberts, M. S. E., & Romani, R. W. 2005, 
ApJ, 627, 904 

Ransom, S. M., Eikenberry, S. S., & Middleditch, 
J. 2002, AJ, 124, 1788 

Roberts, M. 2009, talk presented at Super- 
nova Remnants And Pulsar Wind Nebulae In 
The Chandra Era, July 2009, |http://cxc.| 
harvard . edu/ cdo/ snr09/program . html 



57 



Roberts, M. S. E. 2008, in American Institute of 
Physics Conference Series, Vol. 983, 40 Years 
of Pulsars: Millisecond Pulsars, Magnetars and 
More, ed. C. Bassa, Z. Wang, A. Cumming, & 
V. M. Kaspi, 621-623 

Roberts, M. S. E., Hessels, J. W. T., Ransom, 
S. M., Kaspi, V. M., Freire, P. C. C, Crawford, 
F., & Lorimer, D. R. 2002, ApJ, 577, L19 

Roberts, M. S. E., Romani, R. W., & Kawai, N. 
2001, ApJS, 133, 451 

Saz Parkinson, P. M., et al. 2010, ApJ, 725, 571 

Seidelmann, P. K. 1992, Explanatory Supplement 
to the Astronomical Almanac (University Sci- 
ence Books) 

Smith, D. A., et al. 2008, A&A, 492, 923 

Stairs, I. H. 2003, Living Reviews in Relativity, 6, 
5 

Standish, E. M. 1998, JPL IOM 312.F-98-048 

Story, S. A., Gonthier, P. L., & Harding, A. K. 
2007, ApJ, 671, 713 

Taylor, J. H. 1992, Royal Society of London Philo- 
sophical Transactions Series A, 341, 117 

Taylor, J. H., Fowler, L. A., & McCulloch, P. M. 
1979, Nature, 277, 437 

Theureau, C, Parent, D., Cognard, I., et al. 2010, 
A&A, in press 

Trcpl, L., Hui, C. Y., Cheng, K. S., Takata, J., 
Wang, Y., Liu, Z. Y., & Wang, N. 2010, MN- 
RAS, 405, 1339 

Watters, K. P., Romani, R. W., Weltevrede, P., & 
Johnston, S. 2009, ApJ, 695, 1289 

Weisskopf, M. C, et al. 2006, ApJ, 652, 387 

Yadigaroglu, I., & Romani, R. W. 1995, ApJ, 449, 
211 

Zavlin, V. E., Pavlov, G. C, & Shibanov, Y. A. 
1996, A&A, 315, 141 



This 2-column preprint was prepared with the AAS IATgX 
macros v5.2. 



58 



A. Calculation of Position Offsets 

The light travel time delay r (i.e. Roemer delay) across the solar system from a pulsar at ecliptic coordi- 
nates A (longitude) and /3 (latitude) is: 

t — 500 s cos(/3) cos (6{t) + A) , (Al) 

where 0(t) is the orbital phase of the Earth with respect to the vernal equinox. This is an approximate time 
delay since we are assuming that the Earth's orbit is circular. 

If a pulsar is being timed with incorrect ecliptic coordinates such that there exist position offsets AA and 
A/3, there will be a differential time delay At present in the timing residuals: 

At ~ 500 s [cos(/3 + A/3) cos (6(t) + A + AA) - cos(/3) cos(0(i) + A)] . (A2) 

If the positional offsets are small, such that we can use sinx ~ x, cosx <~ 1, and A/3 A A ~ 0, we can use 
trigonometric identities to get: 

At ~ -500 s [AA cos(^) sin(0(i) + A) + A,3 sin(/3) cos (6(t) + A)] . (A3) 

Comparing the trigonometric identity A sin (0(t) + <f>) = Acos(f>sin9(t) + A sin (j> cos 6(t) to the equation 
for At, we see that: 

Acoscp = -500 s AA cos /3 (A4) 
Asincj) = -500 s A^ sin ^, (A5) 

and therefore: 

A cos <b , . . 

The sinusoid amplitude A and phase <j) come from fits to TOA residuals. One way of doing this is by adding 
a binary model to a Tempo2 fit with parameters Al (A=A1) and TO (with PB=365.2424, the solar year; 
OM=0; and E=0 all held fixed in the fit). Then, 

fmod(T0 - 51623.31250, 365.2424) , , IS 
^ 27r 365^424^ +X + a - 

where f mod is the floating point remainder function, 51623.31250 is the MJD of the vernal equinox in 2000 
and a is a correction for Earth's non-circular orbit that can be measured for a particular point in the sky 
by fitting for TO at several values of simulated position offset in the ecliptic longitude direction only. In that 
case, A/3 = 0, and <f> = 0, allowing us to solve for a. 
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